Spin orbit coupling in bulk ZnO and GaN 
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Using group theory and Kane-like k • p model together with the Lowdining partition method, we 
derive the expressions of spin-orbit coupling of electrons and holes, including the linear-fc Rashba 
term due to the intrinsic structure inversion asymmetry and the cubic-fc Dresselhaus term due to 
the bulk inversion asymmetry in wurtzite semiconductors. The coefficients of the electron and hole 
Dresselhaus terms of ZnO and GaN in wurtzite structure and GaN in zinc-blende structure are 
calculated using the nearest-neighbor sp 3 and sp 3 s* tight-binding models separately. 
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I. INTRODUCTION 

The wide gap semiconductors ZnO, GaN and their al- 
loys with wurtzite or zinc-blende structures have pro- 
voked a lot of interest in the last few years, largely 
due to the good optical properties as well as great po- 
tential in optoelectronics.^^ Very recently, some at- 
tention has also been given to the spin properties of 
these semiconductors £ifLL2vL This is partly because of 
the prediction that ZnO can become ferromagnetic with 
a Curie temperature above room temperature if doped 
with manganese^ Ferromagnetism and magnetoresis- 
tance in Co-ZnO magnetic semiconductors have been also 
reported;- Moreover, long spin relaxation time observed 
in these materials is another promising property for possi- 
ble applications in spintronic devices. Ghosh et al& have 
investigated the electron spin properties in n-type ZnO 
structures and discovered the electron spin relaxation 
time varying from 20 ns to 190 ps when the temperature 
increases from 10 to 280 K. Room temperature electron 
spin relaxation as long as 25 ns has also been measured 
by electron paramagnetic resonance spectroscopy in col- 
loidal n-doped ZnO quantum dots3 In GaN, Beschoten 
et aZ.— reported the electron spin lifetime of 20 ns from 
T = 5 K to room temperature. Recently, hole spin re- 
laxation time about 350 ps at 1.7 K in ZnO epilayer has 
also been reported* 1 ^ 

Spin-orbit coupling (SOC) is the key issue of semicon- 
ductor spintronies i 13 ! 14 Most of the proposed schemes of 
electrical generation, manipulation and detection of elec- 
tron or hole spins rely on the SOC. A thorough under- 
standing of the SOC is therefore very important. In con- 
trast to the zinc-blende semiconductors such as GaAs, 
the existence of hexagonal c axis in wurtzite semiconduc- 
tors leads to an added intrinsic wurtzite structure inver- 
sion asymmetry (WSIA) in addition to the bulk inversion 
asymmetry (BIA)j 15 ' 16 Therefore, the electron spin split- 
tings include both the Dresselhaus effect (cubic in k) and 
Rashba effect (linear in fc)j 17 i 18 The Rashba effect has 
been vigorously discussed using group theory and k ■ p 



arguments by Lew Yan Voon et al.*^ It was pointed out 
that the Rashba SOC coefficient is very small (about 1.1 
meV-A) in ZnOJ^ Majewski and Vogl reported a value 
of 9 meV-A of the Rashba coefficient in GaN based on 
the first principle calculation. 20 Magneto-transport mea- 
surements in GaN heterostructure gave the Rashba co- 
efficient ranging from 5.5 to 10.01 meV-A i 21 i 22 ' 23 i 24 For 
electron Dresselhaus effect in wurtzite semiconductors, 
Wang et al. gave the form of electron Dresselhaus SOC 
H so = le{bk1 — k 2 )(a x k v — a y k x ) recently with j e and b 
the SOC parameters. By fitting a spin degenerate surface 
near the T point, they obtained b = 4.028 and -f e ~ 0.74 
meV-A 3 for AINj 2 ^ However, up to now, there is no report 
on the Dresselhaus coefficients in ZnO and GaN. Addi- 
tionally, there is no investigation on the SOC in valence 
bands. 



In this work, by using group theory and k • p method, 
we first construct the 8x8 Kane model for wurzite struc- 
ture, including the s-p z mixing of the lowest conduction 
band as well as the contributions from the remote bands. 
Then we derive the forms of the SOC for both electron 
and hole by Lowdining partition method. 26 The electron 
and hole Dresselhaus SOC coefficients in wurtzite ZnO 
and GaN are then investigated in detail by using the 
sp 3 nearest-neighbor tight-binding (TB) model, which 
was first put forward for wurtzite semiconductors by 
Kobayashi et al., 27 however, without the SOC effect. We 
incorporate the SOC following the approach by Chadi^ 



This paper is organized as follows: In Sec. II, we 
present the model. We start with the Kane-like k • p 
model for wurtzite semiconductors and derive the expres- 
sions for the SOC for electron and hole perturbatively 
up to third order in Sec. IIA. In Sec. IIB, we briefly in- 
troduce the nearest-neighbor sp 3 TB model for wurtzite 
semiconductors. We present our main results in Sec. Ill 
and conclude in Sec. IV. 
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II. MODEL AND HAMILTONIAN 



where 



In contrast to the form given by Chuang and Chang^ 
in which the s-p z mixing of the lowest conduction band 
and the contributions from the remote bands are missing, 
in this section, we first construct the 8x8 Kane- like k • p 
model for wurtzite semiconductors with the s-p z mixing 
properly included. The contributions from the remote 
bands are also considered perturbatively. Then we de- 
rive the expressions of SOC for both electron and hole. 
Finally, we briefly introduce the sp 3 TB model for the 
wurtzite semiconductors which we use for obtaining the 
SOC coefficients. 



A. Kane-like k ■ p model 

The Schrodinger equation relating the periodic part 
Ut/k(r) of the Bloch function and the energy near the 
band edge has the forr^' ! '' 
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Here V(r) is the periodic potential, H so accounts for the 
spin-orbit interaction, and a% with i = x,y, z are the 
Pauli spin matrices. 

Following Chuang and Chang, the basis functions used 
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near the zone center read |ui 
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Here \S), \X), \Y), and \Z) represent the symmetry of 
the band edge states, and the arrows stand for the spin 
orientation. The z direction corresponds to the c axis of 
the wurtzite crystal. From the Cq v symmetry analysis of 
point group for the wurtzite structure, the k • p Hamilto- 
nian in the basis of \ui) with i = 1, • • • ,8 can be written 
as 



8x8 



K z k 2 

2mo 



E c + 



"3**- + 

iB cv i k z k— 

P2.U , 

Plk z + 





V2iA S2 




iB cv i k z k+ 

E v +A ± + 
A 2 +A'ki 



~^ + + 
B cv3 k z k + 









A 2 

-V2iA s 



V2A 3 



Plk z - 
iB cv2 k^ 

i3,k + 

V2 + 

B cv3 k z k„ 

E v + 
B'kl 







V2A3 








y/2iA a 




E c + 

ft2 fc 2 



iB cv i k z 



Pi L, 

Pik z + 








fc- + 



V2" 
iBcvi k z k- 

E v +A! + 
A 2 +A'kl 



B cv3 k z k— 



-V2iA s 



V2A3 



-E„+Ai 
A 2 



V2^ 










P ± k,~ 

iB c „ 2 fc| 

B„, 3 k z k+ 



iQ h. 



B'kt, 



(4) 



in which /c± = k x ± ik y . The energy parameters are 
defined by (5|ff |5) = E c , (X\H \X) = (Y\H \Y) = 
E v + A 1( (Z|flb|^) = K, (Y\H SZ \X) = iA 2 , and 
(Z\H SX \Y) = -{Z\H ay \X) = iA 3 . The interband 
momentum-matrix elements read Pi = — — (S\p z \Z), 
p 2 = -%(S\P*\X) = -^(S\ Py \Y), 'and Q = 
-JjL<Z|p x |JT) = -^-(Z\p y \Y). It is noted that the pa- 
rarneter Q was omitted in Ref. |29 . The inclusion of Q 
gives an improved description of the valence bands. It 



also determines the hole spin splitting^ 2 - The other pa- 
rameters to', B cv i, B CV 2, B CV 3, A', and B' are related 
with the contributions from the remote bands, which 
were also omitted in Ref. H^. A sz originates from the 
s-p z mixing of the lowest conduction band and has the 
same definition as the energy parameter A3. In addition, 
A sz is closely related to the linear Rashba terms of both 
electron and hole i 19 ' 33 This term was also omitted in Ref. 
I29I . As for the energy parameters Ai, A2 and A3, one 
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can relate them to crystal-field split energy A cr and the 
spin-orbit split-off energy A so by 

Ai = A cr , A 2 = A 3 = l&so- (5) 

Using Lowdining partition method^ one can derive 
the SOC terms for electron and hole, which act as In- 
dependent effective magnetic fields. For the states near T 
point, the electron SOC can be described by [a e fi^(k) + 
7 e f2f (k)] • a with 



a e Of (k) = a e {k yi - 

le n°(k) = le (bk 2 z 
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Here a e Clf (k) and 7 e r2,?(k) are the Rashba and Dressel- 
haus terms, respectively. a e and j e are the corresponding 
SOC coefficients, b is a parameter that can be determined 
by the TB calculation. It should be mentioned that a e 
originates from the s-p z mixing of the lowest conduction 
band, whereas 7 e is closely related with k • p interaction 
with the remote bands having Tq symmetry. 

As we know, for zinc-blende structure, both the heavy 
hole (HH) and light hole (LH) belong to the four di- 
mensional Tsv group presentation whereas the spin split- 
off hole (SOH) has the I^ symmetry 31 ! 34 However, for 
wurtzite structure, the combination of crystal-field and 
SOC energies lead to a three-edge structure at the top 
of the valence band. Two of these three edges are of 
symmetry and the remaining one is of Tg^ symmetry. 
The symmetry of the valence bands are, in the order of 
decreasing energy, T 9v (HH), T 7v (LH), and T 7 > v (crystal- 
field split-off hole) for GaN 3 ^ and T 7v (HH), T 9v (LH), 
and T 7 , v (crystal-field split-off hole) for ZnOJ^ 3 ! This 
anomalous ordering in ZnO results from a negative spin- 
orbit splitting. 38 Hereafter, in the description of the SOC 
of holes of ZnO and GaN, we use the subscript "9" , "7" , 
and "7"' from the symmetry representation to label the 
three kinds of holes above. For the lowest conduction 
band, both ZnO and GaN have T 7c symmetry and we 
use the symbol "e" to denote it. The corresponding SOC 
terms are [o^Hf (k) + 7if2f(k)] • <x, which include both 
the linear Rashba term otitlf-(\t), and the cubic Dressel- 
haus term 7jJl| ? (k), with i = e, 9, 7, 7'. Their expressions 
are given in Table fl] 

The Rashba coefficients with i = e, 7, 7' are ex- 
pressed using the k • p parameters: 
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where E g is the energy gap between the bottom of the 
conduction band and the top of the valence band. As 
for Tg v symmetry, there is no linear-fc splitting, namely 
the Rashba coefficient ag = 0, as pointed out in previous 



literature. 39 ' 40 Moreover, it is noted for the expression 
of a 7 [Eq. ©] that the first term dominates due to the 
large band gap in ZnO and GaN. Therefore, due to the 
negative spin splitting energy in ZnO, both a e and a 7 
are negative. However, they are positive for GaN. The 
situation of a 7 i is reversed. 

As for the Dresselhaus effect, both the SOC coefficient 
7i and the parameter b can be determined by the TB cal- 
culation. While it should be mentioned that 7^ is closely 
related with the remote bands with Tg symmetry. It 
is also noted that the electron Rashba and Dresselhaus 
terms of ZnO and GaN have opposite signs. However, it 
is not the case for and T 7 / v holes. For the rg„ holes, 
only the Dresselhaus term exists. 

For the sake of comparison, SOC terms of zinc-blende 
structures 3 -! are also listed in Table HI For T$ v valence 
band, the SOC form reads 78^8* ' J- Here the small 
modifications from the k • p interactions with the remote 
T3 state and the SOC within the conduction band Fg c 
are omitted. The matrices Ji with i = x,y, z are angular 
momentum matrices for angular momentum J = 3/2. 



B. TB model 

To calculate electron and hole spin splittings, we use 
the sp 3 TB model elaborated by Kobayashi et al.^- 
which has been proven to be an effective approach in 
band structure calculations for wurtzite crystals. 27 ' 41 ' 42 
In this model, the local point symmetry is approximately 
tetrahedral, namley Td symmetry rather than C^ v sym- 
metry. In such an approximation, on-site coupling be- 
tween s and p z orbitals, and the small crystal field split- 
tings between p z orbital and the p XiV orbital are both 
neglected. Thus the model has nine independent param- 
eters: the four on-site matrix elements E(s,a), E(j),a), 
E(s,c), E(p,c) (where s and p refer to the basis states, 
and a and c refer to anion and cation), and five transfer 
matrix elements V^sa), V{spa), V{psa), V{ppTt), and 
V(ppa) [where the orientation of the p orbitals is denoted 
by a and 7r, and the first (second) index refers to the s 
or p state of anion (cation)]. For numerical calculations, 
we use the TB parameters taken from Refs. and EH. 

The SOC is included as outlined by Chadii^ We take 
the spin splitting parameters A Zn =0.335 eV, A o =0.0274 
eV, A Ga =0.174 eV, and A N =0.009 eV^^! 44 It should 
be mentioned that only the cubic term, namely the Dres- 
selhaus spin splitting can be calculated by the TB model 
under Td symmetry approximation. Based on the zero 
spin splitting point calculated by the TB model, one can 
conveniently fit the parameter b in the Dresselhaus SOC 
descriptions listed in Table |U Moreover, the SOC coef- 
ficient 7j can also be obtained from the calculated spin 
splitting energy. 
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TABLE I: Electron and hole SOC terms. a,:f2f (k) and ytOf (k) represent the Rashba and Dresselhaus terms with i = e, 9, 7, 7' 
for wurtzite structure and i = e, 8, 7 for zinc-blende structure, and are the Rashba and Dresselhaus SOC coefficients. 
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III. RESULTS 

Since the linear Rashba term has been investigated in 
Ref. fl9L in this section we mainly discuss the cubic Dres- 
selhaus terms of wurtzite ZnO and GaN, calculated from 
the nearest- neighbor sp 3 TB model. In fact, the electron 
Dresselhaus term becomes the dominant SOC compared 
with the linear Rashba term when the electron concen- 
tration is higher than 10 19 cm -3 for ZnO and 10 20 cm -3 
for GaN. (Note in the literature, the electron concentra- 
tions of ZnO and GaN are in the range of 10 17 cm~ 3 
to 10 21 cm~ 3 4 i 24 i 45 ) In order to compare with the zinc- 
blende structure, the cubic case of GaN is also addressed. 



A. SOC in ZnO 

In Fig. [TJa), we show the cubic electron spin split- 
ting energy around T point for different momentum in 
the planes of k z =0 (red thin curves), 0.1 (blue curves), 
and 0.15 A -1 (green thick curves), respectively. In each 
k z plane, the results along T-K (solid curves) and T- 
M (dashed curves) directions are shown separately. In 
contrast to the spin splitting of GaAs at T point one 
can see from the figure that, for small momentum, the 
splittings show isotropy. Moreover, in the case of k z =0.1 
A -1 , there is a zero splitting point at fcii=0.198 A -1 , 
from which one can determine that b in the Dresselhaus 
term is 3.910. Furthermore, we find from our calculation 
that there is a certain variation of b as shown in the inset 
of Fig. [lja). From the inset one notices that b increases 
monotonically with ku . It also implies that for small mo- 
mentum, b keeps constant around 3.855. Moreover, it is 
independent on the angle of the in-plane wave vector k\\ , 
which can be seen from the fact of the same zero splitting 
points between T-K and T-M directions under the same 
wave vector k z . 

Up till now, there has been no report on the value of 
7 e . In this work, the value of this coefficient is extracted 



from 

7e (k) = A£7/(2|flf (k)|), (11) 

with AE standing for the corresponding Dresselhaus spin 
splitting energy. The results for different momentum are 
shown in Fig. [IJb-d). One can see from the figure that 
for small momentum, 7 e is almost a constant (about 
0.33 eV-A 3 ) with its value being almost two orders of 
magnitude smaller than the electron Dresselhaus coeffi- 
cient 17.0 eV-A 3 in GaAs, 46 thanks to the large band gap 
(E g - 3.3 eV at 300 K) in Zn041 When the momentum 
lies far away from the T point, our calculation indicates 
that the value of j e decreases, either with the increase of 
the in-plane wave vector fcn in Fig. [TJb), or with the out- 
of-plane wave vector k z in Fig. [He). This fc-dependence 
of 7 e is due to the correction of the higher order SOC 
terms for large momentum. Moreover, angle dependence 
of 7 e becomes remarkable for large momentum, which is 
also because of the correction of higher order SOC terms. 

In order to verify the accuracy of our computation, 
we compare our results with the ab initio calculations 
based on atomic sphere approximation linear-muffin-tin- 
orbital (ASA-LMTO) method, as shown in Fig. H In 
the figure, the Dresselhaus splitting is from our compu- 
tation and the Rashba splitting is calculated by using 
the Rashba coefficient extracted from ASA-LMTO cal- 
culations by Lew Yan Voon et al.*£ From Fig. [21 one 
can see that the sum of these two terms is in very good 
agreement with the ASA-LMTO results, which confirms 
the soundness of our Dresselhaus spin splitting calcula- 
tions. This agreement also confirms that the Rashba and 
Dresselhaus coefficients of electrons have opposite signs. 

We now turn to the case of hole SOC. The spin split- 
tings for holes with Tg Vl T-? v , and Tt v symmetries are 
plotted against in-plane momentum in Fig. OJa) and (b) 
in the planes of k z — and 0.04 A" 1 respectively. One 
finds from Fig.OJa) that as fcii < 0.3 A -1 , the spin split- 
tings for the holes with Tj v and Tj' v symmetry increase 
drastically and become much larger than that of the holes 
with Tg v symmetry. The corresponding SOC coefficients 
7i, with i = 9, 7, 7' are also calculated similar to the case 
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FIG. 2: (Color Online) Electron spin splitting in ZnO against 
momentum along T-K direction. The splitting from the 
Rashba term and from the ASA-LMTO calculation are taken 
from Ref. GJ. 



TABLE II: Rashba (at) (in meV-A) and Dresselhaus (ji) (in 
eV-A 3 ) SOC coefficients in wurtzite ZnO and GaN at small 
momentum with i = e, 9, 7, 7 . 
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6.1 
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- 45 b 


32 b 


0.32 
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FIG. 1: (Color Online) (a) Electron spin splitting of ZnO 
against in-plane momentum along T-K (solid curves) and F- 
M (dashed curves) directions. Red (thin) curves: k z —0; Blue 
curves: fc z =0.1 A -1 ; Green (thick) curves: fc z =0.15 A -1 . In- 
set: The coefficient b as function of k\\ . (b) 7 e vs. k\\ along 
T-K direction under different k z . (c) 7 e vs. k z under different 
k\\ . (d) 7e vs. angle of the in-plane momentum at different k\\ 
in the plane of k z =0. 



of 7 e . The results are shown in Fig. (S^c-e). Again these 
coefficients are constants close to the T point, where 77 
and 77/ are close to each other, i.e., 77 = 6.3 eV-A 3 and 
77' =6.1 eV-A 3 . While for the hole with Tg v symme- 
try, the coefficient 79 is only about 0.09 eV-A 3 . More- 
over, 7^ are also angle- independent for small momentum, 
but become markedly angle-dependent for large momen- 
tum. The corresponding results are shown in Fig. [3Jd) 
and (e). Additionally, from the zero spin splitting point 
in Fig. 02Jb) for the holes with and r7'„ symmetries 
(around 0.098 A -1 in the figure), the parameter b can 
also be extracted. We obtain almost the identical value 
as that obtained from the electron SOC under the same 
momentum. 



B. SOC in GaN 

GaN shares the same electronic wurtzite structure with 
ZnO, and also has a direct wide gap about 3.4 eV4£ In 



a from Ref. Q1; b from Ref. 

Fig- Eta), we show the electron Dresselhaus spin splittings 
around T point for different momentums in the planes of 
k z =0 (red thin curves), 0.1 (blue curves), and 0.15 A -1 
(green thick curves), respectively. One can see from the 
figure that the momentum dependence of the Dressel- 
haus splitting is very similar to that in ZnO. In GaN, 
the Rashba coefficient is also small (about 9 meV-A). We 
can also estimate that, the Dresselhaus effect becomes 
dominant when the electron concentration is larger than 
10 20 cm -3 . From the point of zero spin splitting, one 
can calculate b in the Dresselhaus term. The results are 
shown in the inset of Fig.[3]Ja). For small momentum, one 
finds b ~ 3.959 in GaN and 3.855 previously in ZnO, both 
are different but close to 4.028 — the "universal" value for 
all the wurtzite materials predicted by Wang et al. from 
their A1N investigation. 2 ^ Moreover, with the increase of 
the momentum, similar to the case of ZnO, b can not be 
regarded as a constant and shows marked ^-dependence. 
In addition, it should be mentioned that, for large mo- 
mentum, the universality of b for all the wurtzite ma- 
terials does not hold and this difference enhances with 
the increase of momentum. For example, when k\\ = 0.3 
A -1 , the values for ZnO and GaN are 4.039 and 4.326, 
respectively. 

For the hole spin splitting, one also obtains the similar 
results as those of ZnO, as shown in Fig. |DJb) and (c). 
From Fig. 2J it is implied that both the electron and hole 
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FIG. 3: (Color Online) (a-b) Hole spin splitting in ZnO 
against momentum along Y-K (solid curves) and Y-M (dashed 
curves) directions with fc z =0 in (a) and fc z =0.04 A -1 in (b). 
(c) Hole SOC coefficients 7* (i = 9,7,7') vs. ku along Y-K 
direction; (d-e) 7» vs. angle of the momentum in the plane of 
k z =0. fe|| =0.1 A -1 in (d) and fey =0.2 A* 1 in (e). Red curves 
[thick curves in (a) and (b)]: Yg v ; Blue curves [thin curves in 
(a) and (b)]: Yr v ', Green curves: Y 7 i v . 



spin splittings are isotropic at small momentum. The cor- 
responding SOC coefficients 7j with i = e, 9, 7, 7' as func- 
tions of in-plane momentum are shown in Fig. [5Ja) with 
k z = 0. We find 7 e and 79 are about 0.32 and 0.07 eV-A 3 
for small momentum, respectively. These values are very 
close to the case of ZnO, which we believe are due to their 
similar electronic structures and almost equal band gaps. 
However, as for 77 and 77/, they are about 15.3 and 15.0 
eV-A 3 separately, two times larger than those in ZnO. 
Similarly, 77 and 77/ are close to each other. In addi- 
tion, as the momentum increases, the dependence of the 
SOC coefficients on the angle of the in-plane-momentum 
becomes more and more pronounced, as shown in Fig. 
OJb-d). The Rashba (an) and Dresselhaus (7^) SOC co- 
efficients in wurtzite ZnO and GaN for small momentum 
are summarized in Table [TXJ 

In order to have a comparison with the case of zinc- 
blende structures, we also calculate the spin splitting en- 
ergy in cubic GaN using sp 3 s* TB model based on the 
parameters by Reilly et al.^ which has proven its capa- 
bility in the T valley band structure calculations. 46 Both 
the electron and hole SOCs at the T point in zinc-blende 



FIG. 4: (Color Online) (a) Electron spin splitting in wurtzite 
GaN against momentum along Y-K (solid curves) and Y-M 
(dashed curves) directions. Red (thin) curves: k z —0; Blue 
curves: k z —0.1 A -1 ; Green (thick) curves: fc z =0.15 A -1 . In- 
set: The coefficient b as function of k\\. (b-c) Hole spin split- 
ting in wurtzite GaN against momentum along Y-K (solid 
curves) and Y-M (dashed curves) directions with k z =0 in (b) 
and fc 2 =0.04 A^ 1 in (c). Red (thick) curves: Yg v ; Blue (thin) 
curves: Yt v ; Green curves: Y^'v 

crystals have been extensively discussed in Ref. [U and 
it is known that only the Dresselhaus term exists in this 
structure. 

The numerical results are shown in Fig. [6l Two typi- 
cal directions (along Y-K and T-VF) of spin splittings for 
electron and hole are shown in Fig. [f^a). It is implied 
that even though for small momentum, the spin split- 
tings also show anisotropy. This is not in the case of 
wurtzite structure. Moreover, the anisotropic spin split- 
ting is more distinct for the HH, as shown in Fig.[^a). Of 
course, for large momentum, the spin splittings of both 
structures show anisotropy. In addition, it is interest- 
ing to see that although the lowest conduction band and 
the SOH band have different symmetries (rg c and T7„ 
respectively), they share the same form of SOC. From 
the red and black curves in Fig. (Ha), one can see that 
electron and SOH have similar momentum dependence 
of spin splitting. 

The electron SOC coefficient is shown in Fig.^b) and 
(c). Similar to the wurtzite case, j e keeps constant near 
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FIG. 5: (Color Online) (a) The SOC coefficients 7; with i = 
e, 9, 7, 7' vs. fen along Y-K direction with fc z =0; (b-d) 7; vs. 
angle of the momentum in the plane of fc z =0 with k» =0.1 A -1 
(b), fe|| =0.2 A" 1 (c), and k\\=0.5 A" 1 (d). 



the r point and is independent on the direction of the 
in-plane momentum. Here we obtain j e ~ 0.51 eV-A 3 , 
which is also nearly two orders of magnitude smaller 
than that in GaAs4& The big difference comes from the 
smaller band gap of GaAs. When ku > 0.03 A -1 , 7 e 
decreases drastically with the momentum, and its angle 
dependence also becomes very pronounced. 

The corresponding hole SOC coefficients are shown in 
Fig. [S^d) and (e). Near the T point, we find 73 and 77 
around 30 and 25 eV-A 3 , respectively. Additionally, 7§ 
shows stronger angle dependence of the momentum than 
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FIG. 6: (Color Online) (a) Electron and hole spin splitting in 
zinc-blende GaN against momentum along Y-K (solid curves) 
and Y-W (dashed curves) directions. Red curves: electron; 
Blue curves: HH; Green curves: LH; Black curves: SOH. (b) 
The corresponding electron SOC coefficient y e vs. k« along 
Y-K direction, (c) 7 e vs. the angle of the momentum at 
different fen (k z = 0). (d) 78 (red curve) and 77 (blue curve) 
vs. momentum along Y-K direction with fc 2 = 0. (e) 78 (red 
curve) and 77 (blue curve) vs. the angle of the momentum at 
different fey. Solid curves: fey =0.05; Dashed curves: fey =0.1 

A- 1 . 



In conclusion, we have investigated the electron and 
hole SOC up to third order in wurtzite semiconductors. 
Due to the intrinsic structure inversion asymmetry in ad- 
dition to the bulk inversion asymmetry in wurtzite struc- 
ture, both the linear Rashba effect and cubic Dresselhaus 
effect exist. We find the relative signs of the Rashba and 
Dresselhaus terms are opposite for electrons but are the 
same for holes in both ZnO and GaN. For holes with 



symmetry, there is no linear Rashba contribution. 
Due to the negative spin splitting in ZnO, the signs of 
the Rashba coefficients for both electrons and holes are 
opposite with those in GaN. 

The Dresselhaus spin splittings in wurtzite ZnO and 
GaN are calculated by the nearest-neig hbor sp 3 TB 
model. We find that the electron Dresselhaus term 
becomes the dominant SOC compared with the linear 
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Rashba term when the electron concentration is higher 
than 10 19 cm- 3 for ZnO and 10 20 cm" 3 for GaN. While 
for the T-j v and T-?/ v holes, the linear term can not be 
neglected due to the relatively larger hole Rashba SOC 
coefficient. The corresponding SOC coefficients 7 2 ; are 
also calculated, which are momentum-independent very 
near the T point. However, at large momentum, they de- 
crease drastically and show strong angular dependence 
of the in-plane momentum. In addition, we confirm that 
the Dressclhaus parameter b is almost a constant for very 
small momentum and can be regarded approximately as 
an universal value for all the wurtzite materials as pre- 
dicted in Ref. l25l . However, for larger momentum, its 
/c-dependence can not be neglected and the value varies 
strongly for different materials. 

In order to compare with the zinc-blende structure, 
the spin splittings in cubic GaN are also addressed, cal- 
culated using sp 3 s* model. The isotropic spin splitting 
near the T point in wurzite structures does not hold in 
zinc-blende structure. 

We believe that our results are useful for the on-going 
study of spin dynamics of ZnO and GaN. After all, due to 



the small SOC in these wide gap semiconductors, the spin 
relaxation time becomes very long and hence more po- 
tential applications may be figured out using these semi- 
conductors. Finally, we point out that due to the lack of 
the TB parameters for ZnO of zinc-blende structure, we 
cannot calculate the SOC coefficients for this structure. 
More investigations are need. 
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